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^ Abstract 

o 

^js. In this work, a novel approach to explain the survival of sungrazing comets within the 

Roche limit is presented. It is shown that in the case of low tensile strength of the 

cometary nucleus, tidal spUtting of the nucleus can be prevented by the reaction force 

caused by the sublimation of the icy constituents. The survival of Comet C/201 1 W3 

'pH (Lovejoy) within the Roche limit of the Sun is, thus, the result of high tensile strength 

of the nucleus, or the result of the reaction force caused by the strong outgassing of the 






o 



52 icy constituents near the Sun. 
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OO 1. Introduction 

o 

OO 

,__! Shortly before its perihelion passage on 16. Dec. 2011, Comet C/201 1 W3 (Love- 

CO joy) was detected by an Australian amateur astronomer. Quickly it turned out that this 

o 

^vq comet is a member of the Kreutz group, comets that approach the Sun within a few So- 

. . lar radii. An extensive observational program was triggered involving all Solar obser- 

> 

•^ vatories in space. Generally it was expected that Comet C/2011 W3 (Lovejoy) would 

X 

^ not survive its perihelion passage at only 1.2 R^un {Rsun- Solar radius). But Comet 

C/2011 W3 (Lovejoy) re-appeared from behind the Sun and disintegrated 1.6 ± 0.2 
days after perihelion (Sekanina and Chodas 2012[l. 



The Kreutz group is a family of comets on similar highly inclined and high-eccentri- 
city orbits, with semi-major axes of ~ 100 AU and periods of ~ 1000 years. These 
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comets are most probably the products of hierarchical breakups relating back to one 



big cometary nucleus with an estimated diameter of ~ 100 km ( jKreutzj fl 888[ |Sekan 
[mil [20021 [Sekanina and Chodas[ [2004) 1. 

More than 2000 sun-grazing comets have been detected by the the Large Angle and 
Spectrometric Coronagraph (LASCO) onboard the Solar and Heliospheric Observatory 
(SOHO), most of them members of the Kreutz group. Almost daily a new comet is 



found by the sun pointing instruments on various Solar observatories (Knight et al 



2010 1. Generally the observations are not accurate enough to determine their orbital 
parameters, however, the inclination of the orbits are indicative for their membership. 
The vast majority of the sun-grazing group comets are small fragments (meter-sized 
to 10-meter-sized objects) that can only be detected very close to the Sun. Almost all 
are annihilated by the intense Solar radiation and particle bombardment within the 
corona, others hit the Sun directly. Recently, the demise of a comet within the lower 



corona could be observed directly for the first time ( Schrijver et al. 2012 1. A few large 
comets of the Kreutz group were observed from ground. The spectacular Comet Ikeya- 
Seki (C/1965 SI) and the Great Comet of 1882 (C/1882 Rl) are the most prominent 
ones. They probably have split from each other during a perihelion passage about 900 
years ago (Marsden 1967|l. Some comets arrive in pairs. Modeling of their dynamical 



behavior of the Kreutz group shows that these comets can split anywhere along their 
orbits ( [Sekanina and Chodas| [2007[ ). They seem to be distributed along the original 
orbit of their progenitor. 

We assessed the size of the nucleus of Comet C/2011 W3 (Lovejoy) following 



Knight et al. (2010 1, who investigated the light curves of Kxeutz group comets during 
their approach towards the Sun. The brightness of these comets peaks around a helio- 
centric distance of 10 to 12 /?,„„. A comet with a radius of 4 m shows a brightness of 
8 mag at 12 Rsun- The peak brightness of Comet Lovejoy was estimated around -4 mag 
(Karl Battams, NRL, 2012|]] This converts to a radius of ~ 1 km for the nucleus of 
Comet C/201 1 W3 (Lovejoy). At the Comet Lovejoy workshop in Boulder Colorado 
(March 21-22, 2012), J. C. Raymond reported Lyman alpha observations of the hydro- 
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gen (water) production of comet Lovejoy using the UVSC spectrograph onboard the 



SOHO spacecraft (for previous comet observations see e.g. |Ciaravella et aL 2010 1. The 
derived estimate of the nucleus size was 0.4 km in diameter shortly before perihelion. 
The few bright Kreutz group comets (observed from ground) are listed, e.g. in 
Knight et al.| ( |2010 i. Their orbital parameters can be found in the small bodies data 
base of JPl|^ All but one of these comets possess perihelia smaller than the Solar 
Roche limit for self-gravitating solid bodies, 



^roche 



^ ^ Psun 



1/3 



1.8 R, 



(1) 



(calculated for a density of the cometary nucleus of p = 500 kg m , which reflects a 



porous mixture of water ice and dust; [Davidsson et al. 2007 Thomas 2009 1. Here, 
Psun is the mean density of the Sun. The surprising fact that these comets survived their 
Solar passages requires additional forces beyond self-gravity to keep them together. 



Papers dealing with the fates of sun-grazing comets (Huebner 1967 Sekanina 2003 



Iseli et al. 2002 Brown et al. 201 1 1 usually investigate processes that lead to the 



demise of the comets near the sun, but they do not investigate why some survive a 
passage deep inside the Roche limit. One obvious force is due to the tensile strength 
of the cometary material. Estimations of the tensile strength of cometary nuclei are 



ranging from high values, ~ 1 kPa (Rickman 1998 Yelle et al. 2004 Richardson 
et al. 2007|l, to very low values, ~ 1 Pa ( Skorov and Blum 2012|l. Recent observations 



and in particular the Deep Impact results ( A'Heam et al. 2005 Holsapple and Housen 



2007| l have shown that the tensile strength of comets may be very low, clearly below 
1 kPa, most probably in the range of 10 Pa to 100 Pa. Model calculations of the breakup 
of Comet D/1993 F2 (Shoemaker-Levy), when it entered Jupiter's Roche limit, also 
supports the assumption of low cometary tensile strength ( Asphaug and Benz 1994| l. 
Such low values are, however, too small to explain the survival of the brightest Kreutz 
group comets. 

So why do big Kreutz group comets survive the passage within the Roche limit 
but not Comet D/1993 F2 (Shoemaker- Levy)? One obvious diff'erence is that Comet 
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Lovejoy was extremely active near the Sun while Comet D/1993 F2 (Shoemaker-Levy) 
was completely inactive at 5 AU from the Sun. The reaction force caused by the strong 
outgassing (sublimation) of the nucleus near the Sun acts to keep the nucleus together 
and to overcome the tidal disruption. Due to the size dependence of this force (see Eq. 
[8]), it provides an upper limit for the nucleus size. A lower limit is set by the amount of 
eroded material during the perihelion passage. This limit has been widely discussed, 
e.g., by |Huebner| ( [T967) l, |Weissmanl ( [T983) , |Iseh et al.| ( |2002l ) and |Sekaninal ( p003) . The 
actual size of Comet C/201 1 W3 (Lovejoy) probably was close to this lower limit. It 
survived the near Solar passage but its activity terminated in a cloud of dust at 0.4 AU 
(1.6 days after its perihelion, Sekanina and Chodas 20121. 



In this work, we will show a feasibility analysis that sungrazing comets can sur- 
vive a close passage through the Solar corona due to the reaction forces of their high 
outgassing rate. Our model bears many simplifications, so that more detailed analysis 
of the assumptions and effects is necessary in a future study, but provides a sufficient 
criterion for stability, independent of the equation of state of cometary matter. Forces 
acting on sungrazing cometary nuclei are discussed in Sect. l2] In Sect. [3]an upper limit 
for the nucleus size is derived from the balance of forces within the Roche limit as a 
function of the heliocentric distance. Finally, the results are discussed in the concluding 
Sect. H 



2. Forces acting on sungrazing comets 

In this section, relevant forces acting on sungrazing comets are discussed. Without 
better knowledge of Comet C/201 1 W3 (Lovejoy) and for a more general use, we 
assume a spherical, isotropic cometary nucleus at perihelion. Further, we assume that 
all relevant forces can be referred to the midplane of the nucleus (see Fig. [T}. Our 
attempt is to show in a simple model that reaction forces due to the rapid outgassing 
close to the Sun's surface can considerably contribute to the cohesion and, thus, to the 
survival of comet nuclei at close encounters. 

Within the Roche limit, cometary nuclei can be disrupted due to tidal forces acting 
on the material. The tidal force acting on cometary nuclei was investigated in detail by 
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Sun 



Figure 1: Overview of the forces acting on sungrazing comets. 



|Davidsson| ( |200 1 [ ). Here, we adopted his formulation of the tidal force (with = and 
/ = 1 for a spherical nucleus), i.e. 



_ 1 npjM.unR 

r tidal - T ^ 



(2) 

rp 

where Mj,,,, is the solar mass, y is the gravitational constant, p is the mean density of 
the nucleus, R is the radius of the nucleus and rp is the heliocentric distance, or in our 
simplifying case the perihelion distance. 

The self-gravity of the cometary nucleus acts against the tidal forces caused by the 



gravity field of the Sun. Here, the expression derived by Davidsson (2001 with A « 4 



and / = 1 for a spherical nucleus) is used to calculate the self-gravity of the cometary 
nucleus, 

F self -gravity = - 7T p^ y R^ , (3) 

here written as a negative quantity to indicate that this force is trying to prevent the 



tidal splitting. 

The material force further acts against the tidal splitting of the cometary nucleus. 
In our simple case, this force is given by the tensile strength of the material PtensUe 
multiplied by the cross-sectional area of the nucleus at the midplane, 



material 



■^ PtensUe R 



(4) 



Observations and modeling of cometary nuclei yield a broad range of possible values 



for the tensile strength of the material. Richardson et al. (2007i estimated an tensile 
strength of 1 - 10 kPa for the cometary material using results obtained by the Deep Im- 
pact mission. 



Holsapple and Housen ( 2007 1 showed that tensile strengths from Pa to 



1 - 12 kPa could fit the observational data. The tensile strength of the nucleus of Comet 
C/Hyakutake 1996 B2 must be at least 20 Pa to 300 Pa to avoid rotational breakup of the 
nucleus ( |LisseetaL]|1999| l. The breakup of Comet D/1993 F2 (Shoemaker-Levy) yield 
tensile strengths of ~ 270 Pa ( |Greenberg et al.| [T995) and ~ IkPa ( |Rickman||l998| l. 
However, another model calculation of the breakup of Comet D/1993 F2 (Shoemaker- 
Levy) suggest vanishing low tensile strength of the material ( Asphaug and Benz 1994| l. 
Recently, theoretical estimates of the strength of the material composed of individual 
dust and ice aggregates were developed by Skorov and Blum ( 2012| l. They found that 
the tensile strength of hierarchical porous materials can be much lower than that of 
homogeneous (but porous) dust-ice mixtures. For aggregate sizes of 1 mm, they find 
tensile-strength values of only a few Pa. Thus, reasonable values for the tensile strength 
of a microscopically inhomogeneous cometary nucleus composed of dust and water ice 
should fall in the range between ptensUe = ~ 1 Pa (for ~ 1 mm sized aggregates and a 
volume filling factor of (f) - 0.3; |Skorov and Blum| |2012| l and ptensUe = ~ 1 kPa (ho- 
mogeneous dusty bodies at high porosity; Blum et al. 2006| l. 

The sublimation of the icy constituents leads to an additional force acting against 
the separation of the nucleus. This force depends on the amount of absorbed energy 
of the cometary nucleus, given by the energy balance equation for the surface layers 
of the nucleus. We assume that the dominant volatile on the nucleus surface is water 
ice, and that the surface is uniform and is not covered with a substantial dust crust. 
For high surface temperatures {T > 180K; Huebner 1967[ Sekanina and Chodas 



2012|l the absorbed energy is essentially consumed by the sublimation process. The 



heat conduction into the interior of the nucleus and the thermal emission are much 
smaller. In this case, the energy balance equation for the surface of the nucleus is 

^^£™„(rp) = Z(r)A(r). (5) 

The factor 1/4 arises from the assumption that the passing radiative flux is uniformly 
redistributed to the entire surface by the coma around the nucleus. This assumption is 
motivatied by the fact that a coma around a cometary nucleus can significantly influ- 
ence the amount and the distribution of the incident radiation ( |Hellmich]|1981| l. Here, 

E,„„irp)-nFi^t^, (6) 

\ I-P I X2 

is the energy flux density from the Sun at the position of the cometary nucleus, with 
ttF = 6.37xlO^Wm-2,;ri = (1 - [1 - (Ri,J rj)]'/^) /2 andxi = RL/(4rl) 
(Huebner et al. 2007|l. This formulation takes the finite energy flux close to the Sun 



due to the extended solar surface into account. Z(T) and A(r) are the subhmation rate 
and the latent heat of water ice, respectively. The factor ^ was added to the model to 
vary the absorbed energy by the surface layers of the nucleus. The variation of the 
absorbed energy can have different causes: 

1 . The albedo of the surface materials influences the amount of absorbed energy by 
the surface layers. 

2. The presence of an optically thick coma can lead to an increase of the effective 



absorbing cross section of the comet (Hellmich 19811. However, the energy 
distribution inside the coma can either increase, or decrease the radiative flux to 
the surface of the nucleus. 
3. Additional energy input can also be caused by the solar wind, hitting the cometary 
nucleus ( [Brown et al.[|201 1) . 

For water ice, the sublimation rate can be calculated using the Hertz-Knudsen for- 
mula (|Knudsen) [1909}, 



Z{T) ^ a{T)p,„{T)^—-—, (7) 

2nkT 



where a(T) is the subhmation coefficient, m is the mass of a water molecule, k is the 



Boltzmann's constant and PsaiiT) is the saturation pressure of water ice (see Gundlach 
|et al.[|20lT| and references therein for details). 

The resulting force caused by the outgassing of the water ice stems from the mo- 
mentum transfer of the water molecules leaving the surface and can be written as the 
product of the gas pressure and the cross-sectional area of the nucleus at the midplane, 

F outgassing — —^Psat(T)R ■ (8) 

The temperature is calculated from the energy balance given by Eq. E\ 

Simple estimates show that the centrifugal force can be neglected for the slow ro- 
tator and for typical sizes of cometary nuclei. 

3. Size estimation of surviving sungrazing comets within the Roche limit 

With estimates of all forces in place (Eqs. I2]l3]l4j[8]l, the reaction of the comet can 
be estimated both dynamically and statically. A dynamical requirement for survival of 
the comet nucleus can be written in the form 

^ Tidal "•" ^self-gravity "i" ^material ' ^outgassing — '^^cotnet ^esc ^'^ i \^) 

with Meomet, v^sc, c, and T being the comet mass, the escape speed from the comet's 
surface, a factor ~ 1, and the time of perihehon passage, respectively. If Eq. |9]is 
satisfied, then the comet material is unable to leave the comet's Hill sphere during the 
time of passage. A stronger requirement for stability, however, is the condition that the 
comet is even statically stable. This is the case if 

^ tidal "I" ^self-gravity "'" '^material ' '^outgassing — ^- (.^^J 

Here, the cometary material never leaves the comet nucleus surface because the com- 
pressing gravitational, material and outgassing forces exceed the tidal force at any time. 
This criterion is independent of the dynamical constitutive laws of the comet material 
(unlike Eq. l9]l so that it can be applied to comets (whose material properties and con- 
stitutive laws are basically unknown) with only the tensile strength as a free parameter. 



Eq. 10 can be used to derive the maximum radius of a cometary nucleus able to 
sustain the tidal splitting within the Roche limit in the static case. Even for low tensile 
strength of the cometary nucleus, tidal breakup can be prevented by the reaction force 
caused by the sublimation of the icy constituents. Dynamical separation of possible 
fragments due to dynamical evolution is not considered in this approach. 

Fig. |2] shows the resulting maximum radius of the nucleus for different perihelion 
distances (solid curve), using the force caused by the outgassing (Eq. IS] with ^ - 
1) together with the self-gravity force (Eq. l3]l, but neglecting material strength (i.e. 
Fmaieriai - 0). The resulting maximum radius of a cometary nucleus, to survive within 
the Roche limit, due to the material force (Eq. pi and the self-gravity force (Eq. [3]l, 
without the force caused by the outgassing is shown by the dashed-dotted curves, for 

Ptensile - 1 Pa, p,«m7f = lOPn, Ptensile - 1 00 Pa, p,„„,7^ = 1000 Pa (sCC ScCt. |2]l. 

The mean density of the nucleus was assumed to be p = 500 kg m"^, which reflects a 



porous mixture of water ice and dust ([Davidsson et al. 2007 Thomas 2009 |l. 



It is shown in Fig. l2]that the force caused by the outgassing can be orders mag- 
nitude higher than the material force. Thus, active comets are able to withstand tidal 
disruption, even if the tensile strength of the material is negligibly low. 

The lower limit on the radius of the nucleus is given by the thickness of the eroded 
layer during the preperihelion passage and the postperihelion branch of the orbit. The 
thickness of the eroded layer of sungrazing comets was estimated e.g. by |Sekanina| 



( 2003 1 for different perihelion distances and a mean density of the nucleus of p = 
ISOkgm"^. For comparison, these results were recalculated to a mean density of p = 
500 kg m"^ (shown by the asterisks and the dashed curve in Fig. 01. The perihelion 
distance of Comet C/201 1 W3 (Lovejoy) and the Roche limit for a density of p = 
500 kg m"^ are denoted by the dotted vertical lines. 

Using Eq. 10 for Comet C/2011 W3 (Lovejoy; perihelion distance: rp,z,y - 1.27?j„„) 



yields a maximal radius of /?mflA,Ly - 11.0 km (for ^ - I and p,ensiie - Pa). The min- 
imal radius of Comet C/201 1 W3 (Lovejoy) is given by Rmin.u - 0-2 km. A variation 
of the absorbed energy from ^ = 0.5 (50 % less energy absorbed by the surface layers 
of the nucleus) to ^ = 1 .5 (50 % more energy absorbed by the surface layers of the nu- 
cleus) leads to maximum radius of Comet C/201 1 W3 (Lovejoy) of Rmax,Lj - 7.7 km 
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Figure 2: Maximum radius of a cometary nucleus able to survive within the Roche limit as a function of 
the perihelion distance. The solid curve shows the resulting maximum radius using Eq. |10| (with f = 
1), but neglecting the material force. Resulting maximum radii using Eq. |10| and different values for the 
tensile strength, Pfcij,/;, = IPa, p,„,„7p = lOPa, p,„,j;;f, = lOOPa, p,^,,,;/^ = 1000 Pa, but neglecting the 
reaction force caused by the outgassing, are shown by the dasched-dotted curves. For comparison, the lower 
limit on the radius is given by the thickness of the eroded layer during the preperihelion passage and the 
postperihelion branch of the orbit (asterisks and dashed curve; [Sekanina| [2003 } . The dotted vertical lines 
are denoting the perihelion distance of Comet C/201 1 W3 (Lovejoy) and the Roche limit for a density of 
p = 500kgm-3. 



andR„ 



13.6 km, respectively. 



4. Discussion and conclusion 

In this work, we have shown that sungrazing comets are able to survive within the 
Roche limit due the reaction force caused by the outgassing of the icy constituents, 
even if the tensile strength of the material is low {ptensik < 1 kPa). The assumption 
that the entire surface of the nucleus is composed of pure water ice yields an upper 



10 



limit for the radius of the nucleus able to survive within the Roche limit of the Sun, due 
to the self-gravity force and the reaction force caused by the outgassing. 

Using this approach, we derive a maximum radius of R„,ax,Lj - 1 1 -0 km (for ^ - 
i, Piensiie - Pa and p = SOOkgm"-'; see Fig. pb for Comet C/2011 W3 (Lovejoy). 
The variation of the energy available for the sublimation process from ^ - 0.5 to 
^ - 1.5 {see Eq.E\ yields R,„ax,Lj - 7.7km and /?,„(, j^y = 13.6km, respectively. 

The lower limit on the radius of the nucleus is given by the thickness of the eroded 
layer during the preperihelion passage and the postperihelion branch of the orbit. Thus, 
the nucleus of Comet C/2011 W3 (Lovejoy) must have been bigger than Rmin,u - 
0.2 km in order to provide enough material for the erosion. The radius estimated from 
the observed brightness {Ru x 1 km, see Sect. fT) and the radius derived from Lyman 
alpha observations (Ru = 0.2 km, see Sect. fT), are in qualitative agreement with our 
results. 

Due to the outgassing of the icy constituents, the maximum radius of sungrazing 
comets able to survive within the Roche limit is relatively large (see Fig. l2]|. However, 
if the effective gas production decreases, the outgassing force decreases and, there- 
with, the maximum radius of the nucleus able to survive within the Roche limit also 
decreases. Thus, cometary nuclei with low tensile strength can only survive within the 
Roche limit if they are active. Members of the Kreutz group comets like Comet C/201 1 
W3 (Lovejoy) are probably young fragments of a big progenitor comet ( jSekanina and| 
|Chodas| |2002| l and are, thus, active. 

Two very big Kreutz group comets, 1882 II and 1963 V were observed within the 
Roche limit of the Sun (/?ig82// = ~ 31 kmand7?i%3v -~ 14km; Knight etal. 2010|l. 



The perihelion distances were rpj882// - l-67/?5„„ and r^jgesv - 1. 09 7Jj„„, respec- 
tively. Comet 1882 II had broken into at least five fragments (Gill 1883[l. This ob- 



servation is in agreement with our model, because the estimated radius of Comet 1882 
II was bigger than the derived maximum radius for the survival of sungrazing comets 
within the Roche limit (see Fig. l2|. The survival of comet 1963 V can be explained 
with our model within the error of our model and the error of the size estimation. 

A considerably higher activity (caused by, e.g., volatiles with a much higher satu- 
ration pressure; e.g. 20% of the global gas production rate of comet 109P/Hartley was 
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CO2; A'Heam et al. 201 1 1, or additional forces, e. g. the reaction force caused by the 
momentum transfer of the solar wind, can additionally increase the maximum radius 
of sungrazing comets able to withstand tidal disruption within the Roche limit. 

In the case of hemispheric activity, half of the surface receives the total amount 
of the incoming energy. This leads to an increased sublimation pressure and, thus, to 
an increased reaction force caused by the outgassing. However, this resulting force 
only acts on one hemisphere. Thus, the assumption of hemispheric activity does not 
qualitatively change the occurrence of this effect, but can cause a slightly different 
maximum radius of the nucleus. 

Unfortunately, the survival of sungrazing comets within the Roche Umit does not 
provide any information about the tensile strength of the cometary nucleus as long as 
the the reaction force caused by the outgassing exceeds the material force. 

References 

A'Hearn, M.F., Belton, M.J.S., Delamere, W.A., Feaga, L.M., Hampton, D., Kissel, 
J., Klaasen, K.P., McFadden, L.A., Meech, K.J., Melosh, H.J., Schultz, P.H., Sun- 
shine, J.M., Thomas, P.C., Veverka, J., Wellnitz, D.D., Yeomans, D.K., Besse, S., 
Bodewits, D., Bowling, T.J., Carcich, B.T., Collins, S.M., Farnham, T.L., Groussin, 
O., Hermalyn, B., Kelley, M.S., Kelley, M.S., Li, J.Y., Lindler, D.J., Lisse, CM., 
McLaughHn, S.A., Merlin, F, Protopapa, S., Richardson, J.E., WilHams, J.L., 2011. 
EPOXl at Comet Hartley 2. Science 332, 1396-1400. 

A'Hearn, M.F, Belton, M.J.S., Delamere, W.A., Kissel, J., Klaasen, K.P, McFad- 
den, L.A., Meech, K.J., Melosh, H.J., Schultz, PH., Sunshine, J.M., Thomas, PC, 
Veverka, J., Yeomans, D.K., Baca, M.W., Busko, 1., Crockett, C.J., Collins, S.M., 
Desnoyer, M., Eberhardy, C.A., Ernst, CM., Farnham, T.L., Feaga, L., Groussin, 
O., Hampton, D., Ipatov, S.I., Li, J.Y, Lindler, D., Lisse, CM., Mastrodemos, N., 
Owen, W.M., Richai'dson, J.E., Wellnitz, D.D., White, R.L., 2005. Deep Impact: 
Excavating Comet Tempel 1. Science 310, 258-264. 

Asphaug, E., Benz, W., 1994. Density of Comet Shoemaker-Levy 9 deduced by mod- 
elling breakup of the parent 'rubble pile' . Nature 370, 120-124. 



12 



Blum, J., Schrapler, R., Davidsson, B.J.R., Trigo-Rodriguez, J.M., 2006. The physics 
of protoplanetesimal dust agglomerates. I. Mechanical properties and relations to 
primitive bodies in the Solar System. The Astrophysical Journal 652, 1768 -1781. 

Brown, J.C., Potts, H.E., Porter, L.J., le Chat, G., 201 1. Mass loss, destruction and de- 
tection of sun-grazing & -impacting cometary nuclei. Astronomy and Astrophysics 

535, A71. 

Ciaravella, A., Raymond, J.C, Giordano, S., 2010. Ultraviolet spectra of the C-2003K7 
Comet: Evidence for dust sublimation in Si and C lines. The Astrophysical Journal 
713,L69-L73. 

Davidsson, B.J.R., 2001. Tidal splitting and rotational breakup of solid biaxial ellip- 
soids. Icarus 149, 375-383. 

Davidsson, B.J.R., Gutierrez, P.J., Rickman, H., 2007. Nucleus properties of Comet 
9P/Tempel 1 estimated from non-gravitational force modeling. Icarus 191, 547-561. 

Gill, D., 1883. Note on the nucleus of the great comet (b) 1882. Monthly Notices of 
the Royal Astronomical Society 43, 319. 

Greenberg, J.M., Mizutani, H., Yamamoto, T, 1995. A new derivation of the tensile 
strength of cometary nuclei: application to Comet Shoemaker Levy 9. Astronomy 
and Astrophysics 295, L35-L38. 

Gundlach, B., Skorov, Y.V., Blum, J., 2011. Outgassing of icy bodies in the solar 
system - 1. The sublimation of hexagonal water ice through dust layers. Icarus 213, 
710-719. 

Hellmich, R., 1981. The influence of the radiation transfer in cometary dust halos on 
the production rates of gas and dust. Astronomy and Astrophysics 93, 341-346. 

Holsapple, K.A., Housen, K.R., 2007. A crater and its ejecta: An interpretation of 
Deep Impact. Icarus 191, 586-597. 

Huebner, W.F., 1967. Diminution of cometary heads due to perihelion passage. 
Zeitschrift fiir Astrophysik 65, 185. 



13 



Huebner, W.F., Boice, D.C., Schwadron, N.A., 2007. Sungrazing comets as solar 
probes and dust analyzers. Advances in Space Research 39, 413^20. 

Iseli, M., Kiippers, M., Benz, W., Bochsler, P., 2002. Sungrazing comets: Properties of 
nuclei and in situ detectability of cometary ions at 1 AU. Icarus 155, 350-364. 

Knight, M.M., A'Heam, M.F., Biesecker, D.A., Faury, G., Hamilton, D.P., Lamy, P., 
A. Llebaria, M.M., AHearn, M.F, Biesecker, D.A., Faury, G., Hamilton, D.P, 
Lamy, P., Llebaria, A., 2010. Photometric study of the Kreutz comets observed 
by SOHO from 1996 to 2005. The Astronomical Journal 139, 926-949. 

Knudsen, M., 1909. Die Gesetze der Molekularstromung und der inneren Rei- 
bungsstromung der Gase durch Rohren. Annalen der Physik 333, 75-130. 

Kreutz, H.C.F, 1888. Untersuchungen iiber das Cometensystem 1843 I, 1880 I und 
1882 II. I. Theil. Publication der Koeniglichen Sternwarte in Kiel . 

Lisse, CM., Fernandez, Y.R., Kundu, A., A'Heam, M.F., Dayal, A., Deutsch, L.K., 
Fazio, G.G., Hora, J.L., Hoffmann, W.F., 1999. The nucleus of Comet Hyakutake 
(C/1996 B2). Icarus 140, 189-204. 

Marsden, B.G., 1967. The sungrazing comet group. Astronomical Journal 72, 1170- 
1183. 

Richardson, J.E., Melosh, H.J., Lisse, CM., Carcich, B., 2007. A ballistics analysis of 
the Deep Impact ejecta plume: Determining Comet Tempel 1 's gravity, mass, and 
density. Icarus 190, 357-390. 

Rickman, H., 1998. Solar System Ices. Kluwer Academic, The Netherlands, chapter 
Composition and physical properties of comets. 

Schrijver, C.J., Brown, J.C, Battams, K., Saint-Hilaire, P., Liu, W., Hudson, H., Pes- 
nell, W.D., 2012. Destruction of sun-grazing Comet C/2011 N3 (SOHO) within the 
low Solar corona. Science 335, 324-328. 



14 



Sekanina, Z., 2002. Statistical investigation and modeling of sungrazing comets dis- 
covered with the Solar and heliospheric observatory. The Astrophysical Journal 566, 

577-598. 

Sekanina, Z., 2003. Erosion model for the sungrazing comets observed with the Solar 
and heliospheric observatory. The Astrophysical Journal 597, 1237-1265. 

Sekanina, Z., Chodas, P.W., 2002. Fragmentation origin of major sungrazing comets 
C/1970K1,C/1880 CI, and C/1843D1. The Astrophysical Journal 581, 1389-1398. 

Sekanina, Z., Chodas, P.W., 2004. Fragmentation hierarchy of bright sungraz- 
ing comets and the birth and orbital evolution of the Kreutz system. I. Two- 
superfragment model. The Astrophysical Journal 607, 620-639. 

Sekanina, Z., Chodas, P.W., 2007. Fragmentation hierarchy of bright sungrazing 
comets and the birth and orbital evolution of the Kreutz system. II. The case for 
cascading fragmentation. The Astronomical Journal 663, 657-676. 

Sekanina, Z., Chodas, P.W., 2012. Comet C/2011 W3 (Lovejoy): Orbit determination, 
outbursts, disintegration of nucleus, dust-tail morphology, and relationship to new 
cluster of bright sungrazers. The Astrophysical Journal 757, 127. 

Skorov, Y.V., Blum, J., 2012. Dust release and tensile strength of the non-volatile layer 
of cometary nuclei. Icarus 221, 1-11. 

Thomas, N., 2009. The nuclei of jupiter family comets: A critical review of our present 
knowledge. Planetary and Space Science 57, 1106-1117. 

Weissman, PR., 1983. Cometary impacts with the Sun: Physical and dynamical con- 
siderations. Icarus 55, 448^54. 

Yelle, R.V., Soderblom, L.A., Jokipii, J.R., 2004. Formation of jets in Comet 
19P/Borrelly by subsurface geysers. Icarus 167, 30-36. 



15 



